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Abstract 

In relativistic heavy ion collisions, top quarks can be produced by photon- 
gluon fusion when a photon from the Weizsacker- Williams virtual photon field 
of one nucleus interacts with a gluon in the other nucleus. Photoproduction 
with heavy ions at the Large Hadron Collider (LHC) will be the first accessible 
non-hadronic top production channel. We calculate the it photoproduction 
cross sections, pair mass and top quark rapidity distributions in peripheral 
lead- lead and oxygen-oxygen collisions. The cross sections are sensitive to the 
top quark charge and the large-Q 2 gluon distribution in the nucleus. We find 
a cross section of 15 pb in oxygen-oxygen collisions, leading to 210 pairs in a 
one month (10 6 s) LHC run. In pA collisions, the rate is higher, 1100 pairs 
per month for pO. A comparison of the A A and pA data might allow for a 
study of gluon shadowing at high Q 2 . 
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Due to their large charges, relativistic heavy ions carry strong electromagnetic fields 
which may be treated as virtual photon fields. In a relativistic ion collider, these fields inter- 
act with target nuclei in the opposing beam, resulting in high luminosities for photonuclear 
interactions. Because of the large Lorentz boosts, high photon-nucleon center of mass ener- 
gies are reached. Previous studies have considered photoproduction of charm jlj] and bottom 
quarks as well as nuclear breakup || and vector meson production ||. These calculations 
all considered peripheral collisions, with impact parameter, b, greater than twice the nuclear 
radius Ra so that the two nuclei do not interact hadronically. Peripheral collisions will be 
studied experimentally at RHIC || and the LHC |J, now under construction at CERN. 

Here, we consider the production of top quarks via photon-gluon fusion, paralleling 
previous calculations of photoproduction in heavy ion collisions Jl|J^]. We calculate top 
production from the heaviest, Pb, and lightest, O, ions planned for the LHC. In these 
collisions, a ti is produced in the reaction j(k) + A(P) — > t(pi) + tip-i) + X where k is the 
four momentum of the photon emitted from the virtual photon field of the projectile nucleus, 
P is the four momentum of the interacting nucleon in target nucleus A, and p\ and P2 are 
the four momenta of the produced t and t quarks. Note that a phton from the target can 
also interact with a nucleon in the projectile. We work in the center of mass (lab) frame. 
The photons are almost real, with \q 2 \ < (Tic/Ra) 2 - The slight virtuality is neglected. 

On the parton level, the photon-gluon fusion reaction is j(k) + g(x2P) — > t(pi) + i(p 2 ) 
where x 2 is the fraction of the target momentum carried by the gluon. To lowest order (LO), 
the ti production cross section is 
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where the partonic invariants s, ti, and U\ are defined as s = (k + X2P) 2 , £1 = (k — 
Pi) 2 — m 2 = (x 2 P — P2) 2 — m 2 , and Ui = {x 2 P — Pi) 2 — rn 2 = (k — p 2 ) 2 — m 2 with top 
quark mass m = 175 GeV. In this case, s = Ak^x 2 m p where 7 is the Lorentz boost of a 
single beam and m p is the proton mass. Here, et = 2/3 is the expected top quark charge, 
a = e 2 /he is the electromagnetic coupling constant, and a s (Q 2 ) ~ 0.11 is the one- loop strong 



coupling constant evaluated at scale Q 2 = m? + p\ where px is the transverse momentum 
of the produced top quark. The large mass of the top quark prevents toponium production, 
allowing the t and t to be treated as free quarks, even near threshold || . 

The hadronic top production cross section is obtained by integrating Eq. ([!]) over x 2 and 
k, since virtual photons are emitted by the nucleus in a continua of four-momenta, 
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where dN/ dk is the photon flux. The factor of two in Eq. (||]) arises because both nuclei emit 
photons and thus serve as targets. The incoherence of top production guarantees that there 
is no interference between the two production sources [|§. The hadronic invariants can be 
defined for a given photon four momuntum k. If the top quark is detected, the invariants are 
S = (k + P) 2 , T\ = (P — pi) 2 — m 2 , and U\ = (k — pi) 2 —m 2 [T(|. The partonic and hadronic 



invariants are related by s = x 2 S, t\ = Ui, and u\ = x 2 T\. Four-momentum conservation 
at leading order gives x 2min = —U\j (S + T{). In addition to the total cross sections, we also 
present the top quark rapidity distribution and the ti invariant mass distribution. We define 
the t rapidity as y and the t rapidity as y 2 . The top quark rapidity is related to the invariant 
Ti by Ti = —\fSQe~ y . The invariant mass of the pair can be determined if both the t and t 
are detected. The invariant mass is s = M 2 = 2Q 2 (1 + cosh(?/ — ^2))- The minimum photon 
momentum necessary to produce a it pair is /c m ; n = M 2 / {A^x 2 m p ). 

We use the MRST LO gluon distribution xg(x, Q 2 ) |Tl| which is considerably softer than 



the older, flat parameterizations such as xg(x) oc (1— x) n |j2j used in earlier photoproduction 
predictions. Calculations with n = 5 result in cross sections 2—3 times higher than those with 
the MRST LO gluon distribution. Shadowing, the modification of the gluon distribution in 



nuclei, is included via a parameterization based on a fit to data ||13|| . There are two caveats 
regarding shadowing. First, the large Q 2 of top production is outside the upper limit of the 
parameterization, Q^ax = 10 GeV 2 . However, the shadowing parameterization continues 
to evolve beyond this point with changes on the percent level between Q 2 = 10 4 GeV 2 and 



6.4 x 10 GeV . Secondly, the probable impact parameter dependence of the shadowing [14 



has been neglected. 

The photon flux is given by the Weizsacker- Williams formulae. The flux is a function of 
the distance from the nucleus, r, 
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where w = kr/^ and K (w) and K\(w) are modified Bessel functions. The photon flux is 
cut off at an energy determined by the size of the nucleus. In the rest frame of the target 
nucleus, the cutoff is boosted to {2-y 2 — 1)Hc/Ra, or 500 TeV for lead and 1800 TeV for 
oxygen. The tt production cross section, 7p — > tt Jl^] , for a photon with the cutoff energy 



of oxygen is a factor of 7.2 times larger than the cross section at the cutoff energy of lead. 
Thus the difference in the energy cutoffs of lead and oxygen is significant. 

The total photon flux striking the target nucleus is the integral of Eq. (|3]) over the 
transverse area of the target at all impact parameters subject to the constraint that the two 
nuclei do not interact hadronically The numerical result agrees to within 15% of the 

analytical result, 
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the integral of Eq. (|j) for r > 2Ra ~ 2.4A 1 / 3 where w R = 2kB.A/"/- The 15% difference can 
serve as a conservative estimate of the uncertainty on the photon flux; this can be checked 
using other photoproduction and two-photon interactions. 

Although our tt calculation is at leading order, the large mass of the top quark ensures 
faster convergence of the perturbative expansion than for the lighter charm and bottom 
quarks where the next-to-leading order corrections lead to factors of 2 — 3 enhancements 
over the LO cross section. The next-to-leading order tt cross sections are only w 40% larger 



than the LO cross sections in pp interactions ||15|| . We assume the NLO enhancement to be 
similar for photoproduction. 

The expected center of mass energies and average luminosities for O+O and Pb+Pb 
collisions at the LHC are given in Table |I[ The numbers are taken from the latest LHC 



machine studies on LHC ion operation ||16|| . The calculation of the average luminosities 
assumes two experiments, e.g. ALICE and CMS, and a bunch spacing of 125 ns. 

The total cross sections for tt pair production with lead and oxygen beams are 550 pb 
and 15 pb. As Table [I] shows, a 10 6 s (one month) heavy ion run at design luminosity will 
produce 0.2 and 210 pairs respectively. The most significant factor in the rate difference 
comes from the luminosity, ~ 3 x 10 4 times higher for oxygen. 

Shadowing plays little role in the total cross sections since at y — and pt = 0, (X2) ~ 
0.06 for lead and 0.05 for oxygen. These values are in a region where shadowing effects are 
small. In addition, at the large Q 2 required for top production, shadowing effects are further 



reduced by Q 2 evolution [fT3 j . The uncertainty in the large-Q 2 shadowing is important since 



no data are available. However, because of the large luminosity gain, lighter mass ions are 
better for studies of nuclear gluon distributions, reducing any possible effect still further. 

The top quark rapidity distributions, obtained from Eq. (Q), are shown in Fig. 1. The 
calculations assume that the photon is in the field of the nucleus coming from positive 
rapidity so that y < corresponds to k < 7x2^ in the center of mass (lab) frame. If the 
photon is emitted by the target instead of the projectile, the resulting top quark rapidity 
distribution would then be the mirror image of the distribution in Fig. 1 around y = 0. 
(Note that this mirror distribution is equivalent to detecting the t instead of the t. Thus 
the t and t distributions are not symmetric around y — 0.) The total top quark rapidity 
distribution is the sum of the curve in Fig. 1 with its mirror image when both nuclei emit 
photons. Roughly half the production is within \y\ < 1, in the central ALICE acceptance 
and almost all tt production falls within the CMS acceptance, \y\ < 2.4. The tt pair invariant 
mass distributions, da/dM, are shown in Fig. 2. The larger photon energy of oxygen results 
in a broader pair mass distribution. 

The only previous study of it production via photon-gluon fusion in heavy ion collisions 
]I7J used a very different photon flux. The flux was integrated over all r, including the 
nuclear interior, modelling the nucleus as a homogenously charged sphere. This approach 



is flawed because the photon flux inside the nucleus is poorly defined. Ref. [18] shows that 



the photon flux varies by orders of magnitude inside the nucleus, depending on the chosen 
nuclear form factor. More importantly, the calculation is incorrect because the photon 
flux inside the nucleus is much higher than that outside the nucleus. The authors apply 
a correction factor of 1/2 but this factor is inappropriately large. Their photon flux (and 
corresponding cross section) for Pb+Pb collisions is a factor of 10 greater than ours. 

Top quark pairs can be observed via their decays, predominantly it — > W + bW~b, where 
the W decays to iv or qq' . The major background to these channels is likely to be hadropro- 
duction of top in grazing peripheral collisons with slightly smaller impact parameter. The 
two reactions can be separated based on the presence of rapidity gaps in the collision, the 
breakup of the colliding nuclei, and a small multiplicity difference. Because the photon is 
colorless, even if the nucleus breaks up, photoproduction events should also have a rapidity 
gap between the photon-emitting nucleus and the it system. The average multiplicity in pp 
collisions at the LHC is expected to be about 45, slightly higher than the roughly 35 parti- 
cles expected in a 7p collision at typical top-production energies [TJ| . The photoproduction 



multiplicity should be further reduced since much of the photon energy is needed to produce 
the it pair, leaving less energy for produce additional particle production. 

For lead, nuclear breakup can be a complication. A lead nucleus may be excited into a 
giant dipole resonance with probability « 35% • (b/2R A ) 2 . When the resonance decays, the 
nucleus will emit one or more neutrons f20fl . To a good approximation, the excitation and 
the photon-gluon fusion occur independently. Thus the two processes factorize and can be 
studied separately pi] . For oxygen, the excitation probability is very small, so one of the 
interacting ions will almost always remain intact. 

The photoproduction cross sections are large compared with the corresponding pp — > itX 
cross sections. At lowest order, a{pp — > itX) = 41 pb at y/s = 5.5 TeV, 7% of the lead 
photoproduction cross section, while a(pp — > itX) = 80 pb at = 7 TeV, 5.3 times 
the corresponding oxygen photoproduction cross section. Thus, only a moderate hadronic 
rejection factor is required. 

Other backgrounds should be small. Hadronic single or double diffractive production 



without accompanying colored interactions, AA — > AAttX occurs only in a narrow range 
of impact parameters. Diffractive production is also suppressed by the 1/M 4 final-state 
mass dependence. Any single diffractive top production will be at larger rapidities than 
photoproduction which is more central. We have calculated (7(77 — > tt) and found it to 
be negligible. Backgrounds from other photoproduction channels should be significantly 
smaller than at hadron colliders because ai^yp — > QQX)/a( , yp — > X) is much larger than 
a(pp — > QQX) / a{pp — > X) at comparable energies. 

We believe that these criteria should allow for at least statistical separation of events 
containing photoproduction of top in oxygen on oxygen collisions. 

The top photoproduction cross sections are also measurable in pA collisions, as may be 
possible at the LHC [fHIEf. Because the proton and the ion must have the same magnetic 
rigidity and Z = A for the proton, these collisions are at somewhat higher per nucleon 
energies than the corresponding AA collisions. The cost is that the center of mass is no 
longer at rest in the lab. Our pA results are calculated in the equal speed system so that 
the 7 is that of the equal speed system. In this case, the photon flux is calculated using 
the analytical expression in Eq. (^) with wr = (r p + R^jk/'y and r p = 0.6 fm is the proton 
radius. 

The expected nucleon-nucleon center-of-mass energies and luminosities for pA collisions 
are shown in Table ||. No official LHC pA luminosities are yet available |22 |. The values in 



Table [T| are based on estimates in Ref. ||23|| , which were obtained assuming a 125 ns bunch 
spacing and one experiment. 

The pA collisions would allow a measurement of the gluon structure function in free 
protons. A comparison between photoproduction in pA and AA collisions at the same 
energy would provide a straightforward measurement of nuclear gluon shadowing at Q 2 
values far above those currently available. The rates are much higher at the maximum pA 
energies, however, since the cross sections per nucleon are larger than in AA collisions. The 
center-of-mass energy is 40-60% higher, leading to a larger boost in the target rest frame, 
increasing the photon cutoff energy to 1200 TeV for lead, more than a factor of two greater 



than in Pb+Pb, and 3600 TeV for oxygen. In addition, the average photon flux is higher 
in pA because single protons can approach the photon-emitting nucleus more closely than a 
proton in the center of another nucleus. The estimated luminosities for pA interactions lead 
to 40 pairs in pPb and 1100 pairs in pO over a 10 6 s run, as shown in Table [TT]. 

Although the top quarks are produced inside the target nucleus, because of their very high 
boost with respect to the target, they decay well outside the nucleus. For a standard model 



top width of T t =1.5 GeV |24| the top quarks typically travel over 100 fm in the nuclear rest 
frame before decaying. Thus the tt pair acts as a dipole with separation hc/m ~ 10~ 3 fm, 
resulting in a small interaction cross section. Therefore, aside from debris from the target 
nucleon, the target nucleus will be relatively undisturbed, leaving the tt pair with relatively 
few accompanying particles. 

In conclusion, top quark pairs will be produced via photon-gluon fusion in heavy ion 
collisions at the LHC Heavy ion photoproduction will be the first accessible non-hadronic 
production channel for tt pairs. A 10 6 s O+O run will produce 210 tt pairs while pA runs 
will result in higher rates, up to 1100 pairs in pO collisions. These events could at least 
be statistically identifiable on the basis of accompanying rapidity gaps, the presence of an 
intact nucleus, and a slightly smaller multiplicity. The data should allow for measurements 
of the top charge and mass. In addition, the nuclear gluon distribution may be measurable 
at large Q 2 . 

This work was supported in part by the Division of Nuclear Physics of the Office of 
High Energy and Nuclear Physics of the U. S. Department of Energy under Contract No. 
DE-AC-03-76SF00098. 
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TABLES 



Ion Species 


V$NN 


AA Luminosity 


a{AA) 


rtaxe 




(TeV) 


(pb-V 1 ) 


(Pb) 


(per month) 


Oxygen 


7.0 


1.4 x 1(T 5 


15 


210 


Lead 


5.5 


4.2 x 1(T 10 


550 


0.2 


TABLE I. Center of mass energies per nucleon-nucleon collision, \/Snn, 


and luminosities for 


heavy ion collisions at the LHC. 


The tt production cross 


sections and rates in 


a one month (10 6 s) 


AA run are also , 


given. 








Ion Species 


V Snn 


pA Luminosity 


a{pA) 


Rate 




(TeV) 


(pb-V 1 ) 


(Pb) 


(per month) 


Oxygen 


9.9 


3.8 x 1(T 4 


2.9 


1100 


Lead 


8.8 


7.3 x 1(T 7 


55 


40 



TABLE II. Center of mass energies per nucleon-nucleon collision, \/Snn, and luminosities for 
pA interactions at the LHC. The tt production cross sections and rates in a 10 6 s pA run are also 
given. 
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FIG. 1. The top quark rapidity distribution, da/dy, for photoproduced top at the LHC. The 
solid curve is for Pb+Pb collisions while the dashed curve is for O+O collisions. In this plot, the 
photon is emitted by the nucleus with positive rapidity. 
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FIG. 2. The tt pair mass distribution, da/dM, for photoproduced top at the LHC. The solid 
curve shows the results for Pb+Pb collisions while the dashed curve is for O+O collisions. 
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